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We study uj — > ii^-n' decay up to including all orders of the chiral expansion and one- loop 
level of mesons in formlism of chiral constituent quark model. This G-parity forbidden decay is 
caused by rriu 7^ rrid and electromagnetic interaction of mesons. We illustrate that in the formlism 
' both nonresonant contact interaction and p resonance exchange contribute to this process, and the 

I contribution from p resonance exchange is dominant. We obtain that transition matrix element is 

^ ^ . < p|i/ptj|aj >= [—(3956±280) — (1697± 130)i]MeV^ , and isospin breaking parameter is md — mu = 

' 3.9 ± 0.22MeV at energy scale ^ m^. 

^ \ 12.39.-x,12.40.Vv,13.25.Jx,12.15.Ff. 

(N : 

(N . 

^ ■ I. INTRODUCTION 

> : 

' The light current quark masses are baisc input of quantum chromodynamics(QCD). The inequality of the light 
, quark masses, especially, m„ 7^ m^, breaks the isospin symmetry or charge symmetry This breaking of isospin 
' symmetry induces various measurable physics processes such as 7r° — 77, A — YP mixing and u — > tt+tt" decay etc. In 
' this paper, we will focus on the lo — > tt+tt" decay, which is considered as the important source of charge symmetry 
breaking in nuclear physics. 

I In ref. |^(we will quote it as I hereafter) we have shown that the chiral expansion at vector meson energy region 
(-H , converges slowly. Therefore, a well-defined effective field theory describing the physics in this energy region must be 
O 1' available for calculation on high order terms of the chiral expansion and meson loops. Obviously, method of chiral 
' , perturbative theory(ChPT) |^ is impractical to cupature the high order term contribution because the number of 
free parameter increases very rapidly as perturbative order raising. In I, following the spirit of Manohar-Georgi(MG) 
(— I \ model Q, we have constructed chiral constituent quark model(ChCQM) including lowest vector meson resonances. 
• • ' The advatanges of this approach are that high order contribution of the chiral expansion and expansion can be 
. !^ \ calculated consistently, and only fewer free parameters are required. Low energy limit and unitarity of the model 
^ . are also examined successfully. In particular, it is a attractive property that although the leading order theoretical 
$^ ' prediction does not match with experimental data, larger contribution from high order of the chiral expansion and 
5^ 1 pseudoscalar meson one-loop corrects theoretical prediction close to data very much. It is just the characteristic of 
the chiral expansion in this energy region. Therefore, in this the present paper, we also need to calculate uj —t tt+tt^ 
decay to include all high order contribution and pseudoscalar meson one-loop correction. 
This research is also motivated by the following reasons: 

i) In the most recent references |5|-|^], the u tt+tt" decay was treated as being dominant via p-resonance exchange, 
and the direct wtt+tt" coupling is neglected. It has been pointed out in reL Q that the neglect of uj "direct" coupling 
to TT^n^ is not valid. It can be naturally understood since tttt can make up of vector-isovector system, whose quantum 
numbers are same to p meson. Thus in an effective lagrangian based on chiral symmetry, every p field can be replaced 
by TTTT and does not conflict with symmetry. Although authors of ref. |^ also pointed out that the present experimental 
data still can not be used to separate "direct" wtttt coupling from uj — p mixing contribution in model-independent 
way, it is very interesting to perform a theorectical investigation on "direct" wtttt coupling contribution. We will show 
that the contribution from interfernce of "direct" unrir coupling and w — p is about 15%. Thus "direct" wtttt coupling 
can not be neglected indeed. 
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ii) The present study involves the investigation of p" — w mixing, which has been an active subject [5- 14]. The mixing 
amplitude for on-mass-shell vector mesons has been observed directly in the measurement of the pion form-factor in 
the time-like region from the process e+e^ tt+tt^ djll]- For roughly twenty years, p'^ — lo mixing amplitude was 
assumed constant or momentum indepcndt even if p and lu have the space-like momenta, far from the on-shell point. 
Several years ago, this assumption was firstly questioned by Goldman et. al. and the mixing amplitude was 
found to be significantly momentum dependent within a simple quark loop model. Subsequently, various authors 
have argued such momentum dependence of the p^ — u mixing amplitude by using various approaches . In 
particular, the authors of ref. [ |l3| has pointed out that p^ — uj mixing amplitude must vanish at = 0(where 
denotes the four-momentum square of the vector mesons) within a broad class of model. This point will be also 
exmined in ChCQM. 

iii) It has been known that uj — )■ tt+tt^ decay amplitude receive the contribution from two sources: isospin symmetry 
breaking due io u ~ d quark mass difference and electromagnetic interaction. In I we have shown that VMD in 
meson physics is natural consequence of the present formlism instead of input. The vector e^e~ decays are also 
predicted successfully. Therefore, the dynamics of electromagnetic interactions of mesons has been well established, 
and the calculation for u> tt+tt" decay from the transition w — > 7 — > p ^ tttt and "direct" w — > 7 — > tttt is 
straightforward. In this the present paper, we will pay our attention to isospin breaking due to m„ 7^ m^. It is 
another purpose of this paper to determine isospin breaking parameter 5mq = rn^ — m„ via to — s- tt+tt^ decay. This 
parameter is urgently wanted by determination of light quark mass ratios. 

The contents of the paper are organized as follows. In sect. 2 we review the basic notations of the chiral constituent 
quark model with the lowest vector meson resonances. In sect. 3, the tree level effective lagrangian, which including 
all order contribution of the chiral expansion, is obtained. The pseudoscalar meson one-loop corrections are calculated 
in sect. 5. In sect. 6, the formulas and numerical results of w — mixing amplitude and lu — > tt+vt^ are given. The 
sect. 7 is devoted to a brief summary. 



II. CHIRAL CONSTITUENT QUARK MODEL WITH VECTOR MESON 

The simplest version of chiral quark model which was originated by Weinberg and developed by Manohar 
and Georgi Q| provides a QCD-inspired description on the simple constituent quark model. In view of this model, 
in the energy region between the chiral symmetry spontaneously broken (CSSB) scale and the confinement scale 
(Aqcd ^ 0.2 — 0.3GeV), the dynamical field degrees of freedom are constituent quarks (quasi-particle of quarks), 
gluons and Goldstone bosons associated with CSSB(these Goldstone bosons correspond to lowest pseudoscalar octet). 
In this quasiparticle description, the effective coupling between gluon and quarks is small and the important interaction 
is the coupling between quarks and Goldstone bosons. In I we have further included the lowest vector meson resonances 
into this formlism. At chiral limit, this model is parameterized by the following chiral constituent quark lagrangian 
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= iqip + f + g^f 75 - ^y)q - mqq +—< V^[/V^C/T > +-to^ < V^F^ > 



(1) 



Here < ... > denotes trace in SU(3) flavour space, q = (qu^qd^ls) are constituent quark fields. denotes vector 
meson octet and singlet, or more convenience, due to OZI rule, they are combined into a singlet "nonet" matrix 



( L 



V 



(2) 



The and T ^ are defined as follows. 



(3) 



and covariant derivative are defined as follows 



(4) 
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where Ip, = v^^ + a^^ and = Vp, — Ufj, are linear combinations of external vector field and axial-vector field a„, £ 
associates with non-linear realization of spontanoeusly broken global chiral symmetry introduced by Weinberg [[l8[ . 
This realization is obtained by specifying the action of global chiral group G = SU{3)l x 517(3)^ on element ^(<i>) of 
the coset space G/SU{3)^: 

m ^ 9Bm)hH^) = h{<^>)a<^)g[, 9L,9ReG, h{^)eH = SU{3),. (5) 

Explicit form of ^($) is usual taken 

^($) =exp{iA°$°(a;)/2}, [/($)= ^2 (6) 

where the Goldstone boson are treated as pseudoscalar meson octet. The compensating SU{3)^, transformation 
h{^) defined by eq.( ||) is th wanted ingredent for a non-linear realization of G. In practice, we shall be interested in 
transformations of V^, A^, and constituent quark fields under SU{3)^. The q, q transform as matter fields of 
SU(3),, 

q^h{<^)q, q^qh\<^). (7) 

The vector meson fields transform homogeneously under SU(3)^, 

V;^/i($)V;/it($), (8) 

which was suggested by Weinberg and developed further by Callan, Coleman et. al. Since under local G, 

the explicit transformations of external vector and axial-vector fields are 

rp = Vp + a^-^ gii{x)rf,gl^{x) + igR{x)df,gl^{x), (9) 
is SU(3)^, is invariant field gradients and transforms as field connection of SU(3),^ 

A^ ;i($)A^/it($), F^ ^ /i($)F^/it($) + h{<^)d^h\<^). (10) 

Thus the lagrangian( 0) is invariant under Ggiobai x Giocai- 

The several remarks are need here. 1) Note that there is no kinetic term for vector meson fields in C^^. Therefore, 
in this formlism the vector mesons are treated as composited fields of constituent equarks instead of fundamental 
fields. The dynamics of vector meson resonances will be generated via loop effects of constituent quarks. 2) Note that 
there is kinetic term of pseudoscalar mesons in . This is different from some other chiral quark models, in which 
there is no such term. Existing of this kinetic term is consistent with basic assumption of our model, because in this 
energy region, the dynamical field degrees of freedom are both constituent quarks and Goldstone bosons associated 
with CSSB. 3) In C^^ the parameter gA — 0.75 is determined by (3 decay of neutron. It has been pointed out in I 
that this value has included effects of intermediate axial-vector resonances exchanges at low energy. In addition, the 
constituent quark mass parameter m ~ 480MeV has been fitted in I via low energy limit of the model. Such large 
value is requireed by convergence of chiral expansion at vector meson energy scale. 

In this paper, we must go beyond chiral limit for obtaining isospin breaking results. The light current quark matrix 
A4 — diag{TO„, md, ms} can be usually included into external scalar fields, i.e., x — s + ip, where s — Sext + -M, Soxt 
and p are scalar and pseudoscalar external fields respectively. The chiral transformation for x is 

X — ' gBxgi- (11) 

Thus together with ^ and ^\ x and x^ can form SU(3)^, invariant scalar source £,^x^^ and pseudoscalar source 

(^^X^^ ~ S,X^0l5- Then current quark mass dependent lagrangian is written 

- Im^M^ + ix^oq - '^m^M^ - ix^cHq. (12) 

The above lagrangian will return to QCD lagrangian ipMip in absence of pseudoscalar mesons at high energy. So that 
there is a free parameter k which can not be determined by symmetry alonely. From viewpoint of phenomenology, 
this ambiguity is similar to Kaplan-Manohar ambiguity in ChPT It will be studied at elsewhere so that we do 
not further discuss it here. In next section, rigorous calculation will show that our results in this paper is independent 
of K. 



3 



To conclude this section, the ChCQM lagrangian including the lowest vector meson resonances and light current 
quark masses is 

- iq{p + r + <?^^75 - if)q - mqq - ^qi^^M^ + ^x^^Q ~ ^^(^^xC^ - Cx^075<? 

+ ^ < V^UV^'U^ > +^ml < V^V > . (13) 

The effective lagrangian describing interaction of vector meson resonances will be generated via loop effects of con- 
stituent quarks. 

III. LEADING ORDER EFFECTIVE LAGRANGIAN IN N^^ EXPANSION 

In this section, we will derive relevant effective lagrangian via calculating one-loop diagrams of constituent quarks, 
which is at the leading order in expansion. In our calculation, the light current quark masses will be treated 
as perturbation and be expanded to the leading order. Pseudoscalar mesons which are localed in external line 
should satisfy soft pion theorem, i.e., 0(where denotes the four-momentum square of pseudoscalar mesons). 

However, q^ is the four-momentum square of vector mesons, and obviously it is not very small comparing with chiral 
symmetry spantoneously broken scale. Therefore, the higher order terms of q^ in the chiral expansion, have significant 
contribution which can not be neglected. Or in the other words, the chiral expansion at vector meson energy scale 
converge slowly. We will rigorously calculate all g^-dependent term contribution. 

We start with constituent quark lagrangian( |l|), and define vector external source V^{a = 0, 1, • • • , 8), axial- vector 
external source , scalar external source 5*" and pseudoscalar external source as follows 
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(where A"'^, • • • , are SU(3) Gell-Mann matrices and A" = y §)• Then in lagrangian( |^), the terms associating with 
constituent quark fields can be rewritten as follow 

£1 - q{ip - m)q + V^qX^^q + igA^^qX'^ri^q - S'^qX'^q ~ P'^qX'^lr.q. (15) 

The effective action describing meson interaction can be obtained via integrating over degrees of freedom of fermions 

e*-^"" = j VqDqe S '^^'■'^''^''^ =<vac,out\in,vac>v^/^^s,p^ (16) 

where < vac, out\in, vac >y a s p vacuum expectation value in presence external sources. The above path integral 
can be performed explicitly, and heat kernal method |21 2^ has been used to regulate the result. However, this method 
is extremely difficult to compute very high order contributions in practice. In I we have provided an equivalent and 
efficient method to evaluate the effective action via calculating one-loop diagrams of constituent quarks directly. This 
method can capture all high order contributions of the chiral expansion. 
In interaction picture, the equation( FgI) is rewritten as follow 



<o|r,e* |0> 

= > I / d Vl-FTT-:! ■ ■■7:r^^n[Pl,---,Pn)0 (pi - P2 Pn) 

^mi(0) + f z / ^•■•^n„(p„...,p„_,), (17) 

where Tq is time-order product of constituent quark fields, £^ is interaction part of lagrangian( |l5|) , n„(pi, • • • ,p„) is 
one- loop effects of constituent quarks with n external sources, Pi,P2, ■ • ■ ,Pn are four-momentas of n external sources 
respectively and 



4 



n„(pi, • • • ,_p„_l) = J d'^Pn'^nipi,- ■ ■ ,Pn)S'^{Pl - P2 Pn) ■ (18) 

To get rid of all disconnected diagrams, we have 

oo 
n=l 

^i = ni(o), (19) 



(27r)4 {2ttY 



-n„(pi, • • • ,p„^i), (n > 2). 



Hereafter we will call Sn as n-point effective action. 

The 5*1 is tadpole-loop contribution of constituent quarks, which is independent of the purpose this paper. The 
two-point effective action iS'2 has been evaluated in I, 
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52 = ^ y d^x < V.UV^W >+J 
^2 (q) = -l;Aiq^){q'd^, - q^q,) < V^^ {q)V'' {-q) >, (20) 



2 

where kinetic term of pseudoscalar mesons has been renormalized, and A{q^) is defined as follow 

A{q') =9'-^ f dt ■ til t) In (1 - fc|l^). (21) 
Here a universal constant of the model, g, is defined to absorbe logarithmic divergence from quark loop integral 



In I we have fitted g = tt ^ \J Nc/i which satisfy the first KSRF sum rule [2^ rigorously. 



A. Three-point effective action 

Up to the leading order of light current quark masses, there are three kinds of three-point effective action. They are 
made up of by external sources yAA, VVS and l/AP respectively. The effective action with external source yAA 
has been obtained in I, 



S^^^ = -4 / / ^,^''"'B{q^)q^ < V'^(g)[A^(x), A.(x)] >, (23) 



where 



B{q') = -g 
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27r^ Jo Jo 



ii(l-ti)(l-<2)g2 



ln(l-Mi^M^^^)]. (24) 



Let us now calculate three-point effective action with external source VVS. Note that since Vf^ = Vf^ + ieQA^ -t- 
TC^^TT + • • •, for Lu — mixing or wtttt coupling vertices, scalar external source S will reduce to constant matrix. 

d^xd'yd^z ■ V;{x)V^{y)S^ < 0\T{qix)rX^qix)q{y)^^X\{y)q{z)X^q{z)}\0 > 

< Mq)M-q)S > Tr^cMSF{l)YSUl + q)r}, (25) 

where Tr(^^^^ denotes trace taking over color and Lorentz space, Spik) = i{^ — m + ie)^^ is propagator of constituent 
quark fields in momentum space. The direct calculation will give 
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247r2m J J (27r 



where 



6t{l - t) 



Next, we calculate three-point effective action with external source VAP. 

^ (0 Mq)A^{k)P{-q -k)> Tr^,,j^)[SF{l)Yl,SF{l ~ kh^S^il + q)^] 

+ < Mq)P{-q - fc)A,(fc) > Tr^cX)[SF{l)rSF{l " g)75^f + fc)7'^75]}- (28) 

Due to oc dfj^TT + • • • and P oc tt + • • •, for purpose of this paper, the soft pion theroem tell us fc^ — > and 
(k + qY — > 0. In addition, we can find that fcpA^(fc) ^ fc^7r(fc) — > 0. Then performing the loop-integral in eq.( p8| ), 
and employing the above discussion in our calculation, we obtain 

S'^'--^9A J (0i(gl{«l('/^) < mq),A^ik)]Pi-k-,) > 



+a2{q'){q'S^, - q^q,) < [V (q), A'' {k)]P{-k - q) >}, (29) 



where 



^ Jo "'0 

tii2(l-ii)g2 



m2-ti(l-ti)(l-i2)g2^' 
a2(g^) = / ' dti dh , , ^ (30) 



■m2-ti(l-ti)(l-i2)g2" 
is 0{mq) and free pars 

vanish, and if V^j = p|^A^, provide an isospin conservation pim vertex which is order m„ -I- urid and much smaller 

(3) 

that leading order vertex. Thus the contribution from Si^ will be omitted in this paper. 



The third three-point effective action is 0{mq) and free parameter K-dependent. However, if = lo^, 6*3 



B. Four-point effective action 

There is only one four-point effective action relating to cj — )■ tt+tt^ decay. It is made up of by four external source 
yAAS*. As shown in the above subsection, here S reduces to a constant matrix. 

iSi^^ j d^xd'^yd'^zd^w ■ V;{x)Al{y)Al{z)S'^{w) 

X < 0|T{q»7^A"g(a;)q-(y)7-75A^(2/)q-(z)7"75A=q(z)g(u;)A'^q(«;)}|0 > 

2 ( d'^q dH 
~^^J (2^ (27r)4 (27r)4 

x{< V^,{q)A,{k)A,{-k - q)S > Tr(,,L)[SUl)l^SF{l ~ qh^JsSpil - q - A:)7"75] 

+ < F^(q)A,(fc)5A,(-fc-g) >Tr(,.i)[5f(? + 9 + fc)7'' 5^(^-1- fc)7''755|,(07"75] 

+ < V^{q)SA,{k)A,{-k - q) > Tr^,,L)[SF{l + qh^SUl)Yl5SF{l ~ kh'^^s]}- (31) 
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To perform integral of four- moment a I in the above equation and employ indentities in Appendix to simplify result, 
we can obtain 



where 



x{h,{q^) < {W{q), S{x)}A-'{x)A''{x) > +h2{q^) < V (q) A'' {x) S (x) A'^ {x) >}, (32) 

3-2i2i2(l + 2ti)(l-t2)gVw2 



hliq') 



f dti-tl f dt2{l~t2) 

Jo Jo 



[l-tlt2{l~t2)qym^^^ 



h (n^) - i di I dtn f ^ 4(l-tl)[3-4^jt2(l-^2)gVm^] . . 

h2iq)- I dt,-t, / dt2{l-t2) [l_t2^,(l_t,)^2/^2]2 ' (33) 



C. Relevant effective vertices at tree level 



In the following we will give all relevant effective vertice at tree level. Since one-loop correction of pseudoscalar 
mesons will be calculated, we also need to include four-pseudoscalar meson vertices, which will be derived from O(p^) 
effective lagrangian of this formlism. The effective vertices involving electromagnetic interaction have been calculated 
up to one-loop level in I. We will quote them in sect. 5 directly. All effective vertices can be divided into two part: 
one is isospin conservation and author is isospin broken. In addition, we should point out that, so far, all meson fields 
are still non-physical. The physical meson fields can be obtained via the following field rescaling which make kinetic 
terms of pseudoscalar mesons and vector mesons into standard form 

TT^^TT, ^K. (34) 

Since in this paper X-mesons only appear as intermediate states in one-loop diagrams, for sake of convenience, we 
neglect the difference between /tt and (since the results yielded by this difference are twofold suppressed by light 
current quark mass expansion and expansion). 

The — uj mixing vertex, which breaks isospin symmetry, is included in eq. ( [2^ ) 

^^P = ,0^2 I 7^T4e'' ''^o(9')(g'<J^. - q^q.)oj'^{q)p'>''{x) < X^{^M^ + ^^M^) > (35) 

= J;,""^"'^ / ^,e^''-ho{q'){q%.-q,q.)c.^{q)p"''{x) + ---. (36) 

The isospin symmetry unbroken vectors ipcp vertex is included in eqs.( |2^) and ( p3| ) 

4''/;°' - -2/ < K(g)[A^ A''] >, (37) 



where 



b{q') = ^[A{q')+glB{q% (38) 



The isospin symmetry broken vectors ipip vertex is include in eqs.( ^) and ( p2|) 



^(A/=i)__^V^ f^^ .q-x(. a -S a ) 
^V^^ - u^rhngj {2n)^^ ^"^"^ 



x{[ho{q^) + -^glh^{q^)]<{V^{q),iM£,+i'^Me}l^''{x)A%x) > 

-\9Ah2{.q^) < V^{q)A''{x)m^ + ^^M^^)A^x) >}. (39) 
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In particular, define 



we have 



= ^N('Z') + -^Aihliq') ^)], (40) 



= -^^^^^ / ^,e^'^-s{q'm,.-q,q.)^^{q) 

x['K+{x)d''TT-{x) - d''TT+{x)TT-{x)] 

= \ J ^,^''-^b{q'm,. - q,q.)p'^{q) 

X {[K+{x)d''K-{x) - d''K+{x)K-{x)] - [K°{x)d''K°{x) - d" K°{x)K^{x)]} 

^%1<'' = -^^^ / - '^^'^^)^"^('^) (41) 

X {[K+{x)d''K-{x) - d''K+{x)K-{x)] + [K"{x)d'' K"{x) - d^K^ixjK^ix)]} 

^I'^KK^^-y ^,e^''-b{q'){q%.-q,q.)0J^{q) 

X {[K+ixJd^K-ix) - d''K+{x)K-{x)] + [K"{x)d'' K"{x) - d''K°{x)K°{x)]} 

X {[K+{x)d''K~{x) - a''i^+(x)i^-(x)] - [K"{x)d''K°{x) - d^K^ixjK^ix)]} 

Up to O(p^), the tree level four-pseudoscalar meson effective lagrangian has been derived in I. The isospin symmetry 
unbroken four-pseudoscalar vertex is included in the following lagrangian 

(1 - gX) < \/^U\/^'U'^\/^U\/''U'^ >, (42) 



12(47r) 



where we have used = Nc/ (37r^). The isospin symmetry broken four-pseudoscalar vertex is proportional to md — mu, 
which is included in the following lagrangian 

4^'='^ = f So < M{U + C/t) > +^<?i < V^[/V^t/t(X;7t + UM)>. (43) 

I can be found that the eqs.( |35|)-( |4^ ) are free parameter n independent. Moreover, we can see that every 
vector—!- Lpip in eq.( ^l| ) includes an antisymmetry factor {q^5^i, — g^gj^) (where q^ denotes four-momenta of vector 
mesons). Thus the first term of eq.( ^ ) does not contribute to w — *■ ■k~^'k~ decay via pseudoscalar meson loops. This 
antisymmetry factor also constrains that the vertices with one of factors KK, d^Kd^K and Kd^K do not contribute 
to (jj — > TT'^TT^ decay via pseudoscalar meson loops. Then the relevant four-pseudoscalar vertices can explicitly read 
as follows, 

2 2 



4 

72 

J TT 

+ -^{d,,K+d,K- + 9,X°9^i^°) [29^71+9"^^- - (1 - g\)d''Tr-d''n+l 
37rV^ 

^kkZI^ = ^glm{m^K+d^K- - mdK'>d^K'>)i,+ d^n' . (44) 
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IV. ONE-LOOP CORRECTIONS OF PSEUDOSCALAR MESONS 



In this section we calculate one-loop correction of mesons. Because of > >> m^, it can be expected 
that the dominant contribution comes from one-loop digrams of pscudoscalar mesons. In addition, wc can treat pion 
as massless particle but must take ^ 0. This difference is very important, since 7r-loop yields imaginary part 
of T-matrix but K-loop does not at m^j scale. In our calculation, the mass difference between and is also 
neglected. 

There are three kinds of one-loop diagrams correcting to "direct" lutttt couping and uj — mixing(fig. 1 and fig. 
2). 




FIG. 1. One-loop correction to "direct" wtttt couping. 




(a) (b) (c) 



FIG. 2. One-loop correction to a; — r/io" mixing. 

It must be pointed out that, in T-matrices yielded by figure l-(b) and figure 2-(b), the contribution of imaginary 
part is dominant. We have shown in I that it can not ensure unitarity of S-matrix if wc only calculate figure l-(b) 
and figure 2-(b). The unitarity can be ensured through summing over all diagrams in chain approximation(fig. 3 and 
fig- 4). 




FIG. 3. Cliain approximation for "direct" uim couping. 




+ • • • 



FIG. 4. Cliain approximation for u) — p mixing. 

A. Tadpole diagram 

Since pion is treated as massless particle, the nonzero tadpole diagram contribution is yielded by K or r/g mesons(fig. 
l-(a) and fig. 2-(a)). For sake of convenience, here we assume m^g = m^. 

A Correction to ui — p^ mixing 

The tree level up^^p^p is contained in eq.(35), 
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In momentum space, the calculation on fig. 2-(a) is straighforward 



/74 1 - O 
-— — y < A3(A'^A°A1+XA°A° + A"A1A") > . (46) 
(27r)4 B ~rnl+it^ ^ ' ^ ' 



(27r)4 /c2 „ ^2 ^i^A. 

^ ' K ' a=4 

The generators A" of SU(N) obey the completeness relations 

N 



V < A"AA°B > = < AS > +2 < A >< S >, 



a=l 



^ < y'A >< A"B > = 2<AB> <A><B>. (47) 



a=l 



* 16 16, 



Then we have 

V < A3(A°A'^7W + TWA" A" + A"7WA") >= ^ < X^M >= ^(to„ - ma). (48) 
Substituting eq.( ^8| ) into eq.( ^6| ) and performing loop integral, we obtain 

- -t^^^"-^^ J ^,e^'^-M<l')i<l%^ - (49) 

where 

C = 7|^^. A.(±^,./"r(i-£,. (50) 

(47r)^ "^i 2 

Here we define a constant A to absorbe quadratic divergence from meson loop integral. Its value has been determined 
as A = 2/3 in I by OZI rule. 

B Correction to "direct" Cijvrvr mixing 

The isospin symmetry broken uj-Aip couping vertex is included in eq.( p^ ). Expanding eq.( |3^ ) to contain four 
pseudoscalar meson fields, we can obtain 



UTT^f^mJ (27r)4^ -.y^ ,y^^.^a v^,y (27r)4 P - + ie 

8 

X ^{< /2X(a''7rA"[A",a'"7r] + A"[A°,9''7r]9'^7r) > 

a=4 

+ i < hd'^Trd'^TriX'^rM + TWA" A" + A"7WA") >} 

= y^l2^^^ y (2^^ ^^'^ ^('^ '^^ " 

x[7r+(a;)a'^7r^(a;) -a''7r+(x)7r-(x)], (51) 

where /2 = diag{l, 1,0}, s(g^), to see eq.( ^ and eq.( ^) has been used. 
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B. /sT-loop contribution 



Here K-loop denotes that one-loop diagrams in fig.l-(c) and fig.2-(c). In this subsection, T will denote time-order 
product of ii'-meson field. 

A Correction to a; — p*^ mixing 

The u — effective action yielded by K-loop is follow 



lg{K-loop) 



Nc mu-md f d q d I ^ ^ ^ 2 r ^ 

" 3^ m J (2^ (2^ ^^^^ ~ ^'^^'^^^^ "'^ ~ ^"^''^ 

xa;'^((?)pO"(-fi)(/ + qYlP^K{l)^K{l + q), (52) 

where Ak{1) = — + ze)~^ is propagator of if-meson. Integrating over in the above equation and defining 
we have 

= ^^^^^^ / ^6(9^)5(9^)^^(9^)9^(9^^ - 9.9.)a;''(9)A(-9). (54) 
The corresponding effective lagrangian read 

= ^^^^^ / ^e--%^)^(5^)S^(9^)9^(9^V - 9.9.)a;''(9)p°''(a=). (55) 
B Correction to "direct" wtttt couping 

The "direct" lotttt couping effective action yielded by if-loop can be evaluated as follow 

^siK-loo,, ^4^^4^ < o|T{4^r^(.)4t^K^(y) + 4^^;^)(x)4tK^(y)}|o > 

_ jV, m„ - mrf /■ ^^^Mg^2 -.-2^. 2^ _ 2 2 _ (3-ffl)iV 

x(9'<5^. - 9M9>^)(fc • 0^"^n9)^+(-9 - k)^-{k)AK{l)AK{l + 9) 
A^c niu-md^ f d'^q d'^k ro.-.2 j-2l/_2\ 1 ./..2\/i , 



7r2/2 m 

x(9'V - 9^9.)fc''a;''(9)7r+(-9 - k)n-{k), (56) 
where we have taken soft pion limit and 3 — gf^ ~ 3 due to g\ <C 0.3. The corresponding effective lagrangian reads 

"-27r2/2 m % ) - -.(9 )(1 + ^)]Ek(9 ) 

x(9'V - 9^9.)'^''(9)[7r+(a;)9^7r-(a;) - a''7r+(a;)7r-(a;)]. (57) 



C. Chain contribution of vr-loop 

Finally, we calculate chain approximation corrections of 7r-loop in fig. 3 and fig. 4. 
A Correction to "direct" wtttt coupling 
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The effective action yielded by 7r-loop in fig. l-(b) is evaluated as follow 

= -J d^xd^ < o\T{£^^^{x)C4Ay)}\o > 



4:Nc niu-md f d'^q d^k dH , q^N, 



x{l ■ k)ruj^'{q)Ti+{-q - k)TT-{k)A^{l)A^{l + k), (58) 

where we have employed soft pion limit and g\ ^ 3, and A7r(0 = i(P + ie)^^ is propagator of pion. Integrating over 
I in the above equation and defining 



(4^) 

we have 



^Aq') -j£^i^ + fdt- t{l - t) In + '-Argi-mq' - Ami)}, 



(59) 



oU^loor)_ 2jVe mu-mq r d^q d^k ^ q^N, , 

x{q^S^,-q^q,)k-'u^{q)7r+{-q-k)7r-{k). (60) 

In eq.( p9| ), Arg{—1) = — tt has been fitted in I due to requirement of unitarity. The corresponding effective lagrangian 
reads 

3vr2/2 m 7 (27r)4 47r2/^ 

x(gV-9pg.)^"(g)k+(:r)9'^^"(-^) -5^^+(a;)^"(a;)]. (61) 
Comparing eq.( |6l] ) and tree level vertex( 41), we can find that every one-loop in fig. 3 contributes a factor 



-S(g2) = -4/-2(l + g|)E.(,2). (62) 
Thus summing over all diagrams in fig. 3, we obtain 

C = 1^^^^^^ / (04 ^'"■'iT^^'^''^-'' - q,q.)u^{q)[.H^)d'''^-{x) - d'^ ^+ {x)^- {x)]. (63) 

B Correction to uj — p'^ mixing 

If in fig. 2-(b), tree level vertex C^^-^t^ is replaced by £*jr7r which contains all diagram contribution in fig. 3, then 
summing tree diagram and fig. 2-(b) is just chain approximation correction to uj — mixing. The effective action 
yielded by 7r-loop in fig. 2-(b) is evaluated as follow 



1 d^xdS < 0\T{C:^,{x)Cp.,iy)}\0 > 



■ q{T!-loop) 



67r2 m J (27r)4 ' '^'^ ' ^'-^ ' 1 + S(g2) 
Thus chain approximation of fig. 4 yield effective lagrangian as follow 



" ' Kq )s{q )^^(q ) — , , ,r — i^^[q)p {-q}- (64) 



= ^^^^V^ / (0^^''''^^^'^°('^') - (i + 5(g2)) ^(g')^(g')^-(g')>(g'^^- - q.q^)^'iq')p°''i^)- (65) 
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Y. w- (P MIXING AND lo tt+tt" DECAY 

Due to VMD, the eleectromagnetic interaction contributes to lo ^ tt+tt^ decay through — * 7 — )■ p*^ ^ tttt and 
w ^ 7 ^ TTTT. In I we have evaluated ptttt vertex, — 7 mixing vertex and "direct" ^tttt vertex up to one-loop level. 
The "direct" 77r7r vertex reads 



= / ^e^''--FAq')A,{q)[^+{x)d^^'K-{x) 5^7r+(a:)^-(x)], (66) 



where is photon field, F^i^q^) is nonresonant background part of pion form factor. Explicitly, FT^{q'^) reads 

1 + S(g2 



W) = l + ,^^, (67) 



where 



b{q^) ^,2, C(g2)Eo(g^) 



2(1 + 30 l + llC/3 ' 

1 

2/^ 



2 



So(g') = Tl[2S.(g2)-I]^(q2)], (68) 
S('Z^) = [l + l^?^]So(a, 



D{p 



2\ 



l + llC/3' 

——2 {A + / dx-xil~ x) In 1 - ^ ^2 

-^*7r^(p'-4m2)}. 



The complete ptttt vertex reads 



pTTTT 

with 



- / -^e^'-^gpUq'KqX^ - q,q.)p'^{q)[7r+{x)d'^7T~{x) - d'^n+ {x)n~ (x)], (69) 



(1 + 2C)(1 + S(g2))- 



ffp^^C'^') - n , o.Nn :w.2^N - (70) 



Moreover, the complete p — 7 mixing vertex reads 

'^^^ = / (04e^'^"^P7(9')('z''5^. - g^g.)p"ng)A^(a;), (71) 

where 

b (q^)~^^-fV)^\l + ^^^] (72) 

5(1 + ^" ^^^ + 2C^ l + S(g2)J- ^'^> 
Thus due to VMD, the complete uj — j mixing vertex can be obtained directly 

6"^ (27r) 

Eqs.( ^1]) and ( ^3|) will lead to w — mixing at the order of ao.m. through the transition process w ^ 7 ^ p", which 
is 



= -^e J AA_e^^-b,,iq'){qH,. - q^q.)p^^{q)A'^ {x). (73) 



C'ZT = l^e' / (04e''^'^P7(9')('?''^P'' - q,qu)u^''{q)p'''{^)- (74) 
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In addition, eqs.(|62) and ( 73) also lead to "direct" ujittt couping at the order of ac.m. through the transition process 
uj "f TTTT, which is 

xw^(q)[7r+(a;)5'^7r-(a;) - d" tt+ {x)tt- {x)]. (75) 
Eq.( |7|) together with eqs.( |9|), (|5|) and Q give the complete uj — mixing vertex as follow 

1^1, - / (04e^'^-e.p(g2)(g2^^^^ - q,,q,)u,^{q)p^-{x)^ (76) 

where vector meson fields have been normalized to physical fields, and 

e.p(9^) = i%!^^Ii^{5-2/.o(<Z^)(l - ^C) + q%{q^)^{q'WK{q^) T^4^]} + Z^lM^ (77) 



The complete "direct" uj-k-k vertex can be obtained via summing eqs.(^), ( p3| ) and ([TS]), 



C-l.-. = / ^e^'^-"5-..(<z')('7''5p. - qpq.)uj\q){T:^(x)d^T:-(x) - a^^+(x)^- (x)]. (78) 



where all meson fields have been normalized to physical fields, and gunmiq ) is defined as follow 



2^ Nc ruu - md r , 9. f 1 10 



6/-^S;,(g^)[8mV-^6(g^) - ^(1 + ^)]} + ^^.(g^)6p.(9^). (79) 



Thus G-parity forbidden cj —> tt+tt includes a nonresonant background contribution, eq.( |78| ), and p resonance 
exchange contribution(eqs.( 69) and (|76|)). The decay width on u mass-shell is 

i(a;^7r tt ) = t^I 2 TV- r 9u^^7r[mJ\ [I —) ' . (80) 

Using the experimental data B{ll) tt+tt^) = (2.21 ± 0.30)% together with eq.( |8^ ), we obtain 

m„ -mrf = -(3.9±0.22)MeV (81) 

at energy scale p ~ m^^. Here the error bar is from the uncertainty in branch ratio of the process tu tt^tt^. In the 
standard way, the lo — p^ mixing amplitude is 



^^;^n„p(g^) =< J d'^xC^p{x)\p > =^ Tiupiq'^) = q^Qu,p{q^ 



(82) 



The off-shell bj — p'^ mixing amplitude is obviously momentum dependent, and vanished at =0. This is consistent 
with the arguement by O'Connell at. ai. in ref. |]l3| that this mixing amplitude must vanish at the transition from 
time-like to space-like four momentum within a broad class of models. In addition, the value of isospin broken 
parameter ( |8l[) leads on lo mass-shell lu — p^ mixing amplitude as follow 

Ren^p(m2 ) = -(3956 ± 280)MeV^ Imn^p(m5) = -(1697 ± 130)MeV^ (83) 

In ref. j^, the on-shcU mixing amplitude has extracted from the e^e^ — > tt+tt^ experimental data in a model- 
dependent way. In eq.( ^3|), the real part of on-shell mixing amplitude agree with result of ref. ||^. The imaginary 
part, however, is much larger than one in ref. which is around — 300MeV^. It must be pointed out that, in ref. 

the author's analysis bases on a model without "direct" lotttt coupling. Therefore, it is insignificant to compare 
the value of on-shell mixing amplitude of this the present paper with one of ref. |^ . Fortunately, the ratio between 
Lj — !■ TTTT decay amplitude and p — > tttt decay amplitude should be model-indenpendent. This value can test whether 
a model is right or not. The on-shell mixing amplitude in ref. [p| yields 
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R^JP = — = -(0.0060 ± 0.0009) + (0.0322 ± 0.0050)i. (84) 
< 7r+7r~|p > 

The present paper predicts 

^ 'T^g,9^p(mg,) gu^.A^l) ^ _(Q ^ Q 0007) ^ ^0 0331 ± 0.0021)i. (85) 

- mf, + impTp gp^r, (m^ ) 

We can see that this theoretical prediction agree with experimental excellently. 
Moreover, the follows have also been revealed in our studies of this paper: 

i) If we take gui-Kv{(f') = in eq.( |o|), we have B{uj tt+tt^) = (2.56 ± 0.34)%. So that the contribution from 
interference between "direct" ujttt: coupling and uj — mixing is about 15%. The dominant contribution are from 
p-resonance exchange. This conclusion indicates all pervious studies which without "direct" wtttt coupling are good 
approximation even though this neglect is an ad hoc assumption. However, in mechanism of w ^ tt+tt" with "direct" 
omegamr coupling, larger imagnary part of on-shell to — mixing amplitude is allowed, but it is not allowed in the 
mechnism without the direct couping. 

ii) If we do not consider the contributions from one-loop diagrams of pscudoscalar mesons, i.e. setting (g^) = 
^7r(9^) = 0, we obtain B{uj tt+tt") = (2.86 ± 0.47)%. Thus the contribution from one-loop of pseudoscalar mesons 
is about 30% and can not be omitted. This conclusion is consistent with I. In addition, in this case, the on-shell w — p° 
mixing amplitude is about — 4700MeV^. So that we can see that the larger imagnary part of on-shell uj — p^ mixing 
amplitude is yielded by pseudoscalar meson loops. In I, we have shown that this larger imagnary part is required by 
the unitarity of this effective field theory. 

VI. SUMMARY 

In this paper, we study G-parity forbidden lo — > tt+tt^ decay up to one-loop level of mesons. This process is yielded 
by isospin symmetry breaking due to m„ ^ and electromagnetic interaction of mesons. The decay amplitude 
contains two parts of contributions which are from "direct" unrir couping and uj — p^ mixing respectively. In the 
previous studies, the "direct" cjtttt couping is neglected. We show that the "direct" wtttt couping and its interference 
with — p'^ mixing contribute to on-shell decay amplitude about 15% only. It also interprets why the previous studies 
are good approximations even without "direct" wtttt couping. We suggest that the decay amplitude ratio i?^^ should 
be model-independent, and our prediction agree with experimental data excellently. 

The formula of a; — p*^ mixing amplitude is also obtained. Since our calculation is beyond the chiral expan- 
sion (including all orders contribution of the chiral expansion) and one-loop contribution of pseudoscalar mesons is 
considered, the momentum-dependence of the off-shell mixing amplitude is very complicated. However, the mixing 
amplitude also vanishes at = 0. For case of on lo mass-shell, the mixing amplitude emerges larger imagnary which 
is from one-loop contribution of pion and is required by unitarity of this effective field theory. 

In our calculation, all vertices are expanded to the leading order light current quark masses. At this order, the 
decay amplitude yielded by isospin broken is proportional to — to„. The theorectical prediction of isospin breaking 
parameter is m^ — m„ = (3.9±0.22)MeV at energy scale p ~ m^^. This value is important for determining light quark 
masses at vector meson energy scale. 



Appendix 

Here we provide some identities which are used in calculation on four-point effective action. In sect. 3.2 we have 
used (7, k and — fc — g to denote four-momentum square of external source V"^, Aj^ and Ao-. For the purpose of this 
paper, Ai,(A:) reduces to k„n{k). So that due to soft pion theorem we have fc^ ^0, {q + k)^ — > and fc^A^(A:) — > 0. 
Moreover, due to space-like condition of vector meson fields, q^Vp{q) = 0, we have 

{Sf^aq. + Sp,q„) < {y^(g), ^}A''(fc)A^(-fc - g) > = 0, 
{Spaq. + Sp,q,) < V''{q)A''{k)SA^~k - q) > = 0. 

quQakp < {V^{q),S}A''ik)A^-k-q) > ^ -^q.Sp, < {V^iq), S}A'' ik)A^-k - q) >, 
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q.qak^ < V^{q)A''{k)SA^-k ~ q) > ^ -y9.<5^. < ^^(g)A'^(fc)5A"(-fc - q) >, (86) 

k^S^a < {V^{q), 5}A'^(fc)A'^(-fc -q)>^ -qj,,, < {V' {q), SjA'' {k)A^ {-k - q) >, 
k^^S,, < V^iq)A''{k)SA''{-k -q)>^ -qj^. < V''{q)A''{k)SA''{-k - q) >, 

The following integral identities are also used in our calculation 

• - '\l^-.n^+llil-.)q^V ~ '\r^-J+~.yt-^)q^]^ 

ACKNOWLEDGMENTS 

This work is partially supported by NSF of China through C. N. Yang and the Grant LWTZ-1298 of Chinese Academy 
of Science. 



[1] GA. Miller, B.M.K. Nefke ns and I. Slaus, Phys. Rep. 194 (1990) 1. 
[2] X.J. Wang and M.L. Yan, [hep-ph/000115C| , submit to Phys. Rev. D. 
[3] J.Gasser and H.Leutwyler, Ann. Phys. 158(1984) 142; Nucl. Phys. B250(1985) 465. 

[4] A. Manohar and H. Georgi, Nucl. Phys. B234 (1984) 189; H. Georgi, Weak Interactions and Modern Particle Theory 

(Benjamin/Cimmings, Menio Park, CA, 1984) sect. 6. 
[5] Coon,S.A., Barrett,R.C.: Phys. Rev. C36 (1987) 2189. 

[6] H.B. O'Connell, A.W. Thomas and A.G. Wilhams, Nucl. phys. A623 (1997) 559. 
[7] A. Bernicha, G. Lopez Castro and J. Pestieau, Phys. Rev. D50 (1994) 4454. 
[8] S. Gardner and H.B. O'Connell, Phys. Rev. D57 (1998) 2716. 
[9] K. Mahman, H.B. O'Connell and A.G. Wilhams, Phys. lett. B376 (1996) 19. 
[10] R. Urech, Phys. Lett. B355 (1995) 308. 

[11] T. Goldman, J.A. Henderson and A.W. Thomas, Few Body Systems 12 (1992) 123. 

[12] J. Piekarewicz and A.G. Williams, Phys. Rev. C47 (1993) 2462; G. Krein, A.W. Thomas and A.G. Williams, Phys. Lett. 
B317 (1993) 293; T. Hatsuda, E.M. Henley T. Meissner and G. Krein, Phys. Rev. C49 (1994) 452; K.L. MitcheU, P.C. 
Tandy, CD. Roberts and R.T. CahiU, Phys. Lett. B335 (1994) 282; R. Friedrich and H. Reinhardt, Nucl. Phys. A594 
(1995) 406; H.B. O'Connell, B.C. Pearce, A.W. Thomas and A.G. Williams, Phys. Lett. B354 (1995) 14; A.N. Mitra and 
K.C. Yang, Phys. Rev. C51 (1995) 3404; K. Maltman, Phys. Lett. B362 (1995) 11; T.D. Cohen and G.A. Miller, Phys. 
Rev. C52 (1995) 3428; K. Mahman, Phys. Rev. D53, (1996)2563, 2573. 

[13] H.B. O'ConneU, B.C. Pearce, A.W. Thomas and A.G. Wilhams, Phys. Lett. B336 (1994) 1. 

[14] M. Benayoun, S. Eidelman, K. Mahman, H.B. O'ConneU, B. Shwartz and A.G. Wilhams, Eur. Phys. J. C2 (1998) 269. 
[15] L.M. Barkov et. ah, Nucl. Phys. B256 (1985) 365. 

[16] J.J. Sakurai, Currents and Mesons, Univertity of Chicago Press (1969). 
[17] S. Weinberg, Physica A96 (1979) 327. 
[18] S. Weinberg, Phys. Rev. 166 (1968) 1568. 

[19] S. Coleman, J. Wess and B. Zumino, Phys. Rev. 177(1969) 2239;C.G. Callan, S. Coleman, J. Wess and B. Zumino, ibid. 
2247. 

[20] D.B. Kaplan and A.V. Manohar, Phys. Rev. Lett. 56 (1986) 2004. 

[21] J. Schwinger, Phys. Rev. 82 (1951) 664;J. Schwinger, Phys. Rev. 93 (1954) 615. 

[22] R.D. BaU, Phys. Rep. 182 (1989) 1. 

[23] K. Kawarabayasi and M. Suzuki, Phys. Rev. Lett. 16 (1966) 255; Riazuddin and Fayyazuddin, Phys. Rev. 147 (1966) 
1507. 

[24] Particle Data Group, C. Caso et. ah, Europ. Phys. J. C (1998). 



16 



